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ABSTRACT 

Using VLT/FLAMES optical integral field unit observations, we present the first spa- 
tially resolved spectroscopic study of the well known blue compact galaxy Haro 11, 
thought to be a local analogue to high redshift Lyman Break Galaxies. Haro 11 dis- 
plays complex emission line profiles, consisting of narrow (FWHM < 200 kms -1 ) and 
broad (FWHM ~ 200-300 kms -1 ) components. We identify three distinct emission 
knots kinematically connected to one another. A chemodynamical analysis is pre- 
sented, revealing that spatially resolved ionic and elemental abundances do not agree 
with those derived from integrated spectra across the galaxy. We conclude that this 
is almost certainly due to the surface-brightness-weighting of electron temperature in 
integrated spectra, leading to higher derived abundances. We find that the eastern 
knot has a low gas density, but a higher temperature (by ~4,000 K) and consequently 
an oxygen abundance ~0.4 dex lower than the neighbouring regions. A region of en- 
hanced N/O is found specifically in Knot C, confirming previous studies that found 
anomalously high N/O ratios in this system. Maps of the Wolf-Rayet feature at 4686 A 
reveal large WR populations (^900-1500 stars) in Knots A and B. The lack of WR 
stars in Knot C, combined with an age of ~7.4 Myr suggest that a recently completed 
WR phase may be responsible for the observed N/O excess. Conversely, the absence of 
N-enriched gas and strong WR emission in Knots A and B suggests we are observing 
these regions at an epoch where stellar ejecta has yet to cool and mix with the ISM. 

Key words: galaxies: abundances - galaxies: individual (Haro 11) - galaxies : dwarf - 
galaxies: kinematics and dynamics - galaxies:interactions - stars: Wolf-Rayet 



1 INTRODUCTION 

Blue Compact Galaxies (BCGs), due to their high star- 
formation rates and low metallicities, are thought to be 
the nearest analogues to young, starbursting galaxies at 
high redshift. However, unlike their distant counterparts, the 
proximity of BCGs provides us with an opportunity to study 
galaxy evolutionary processes in great detail. 

Haro 11 (ESO 350-IG 038) is a well-known member of 
this galaxy class, having been studied in detail in X-ray, 
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sive BCG that consists of three main knots (Fig. [TJ, each 
containing numerous super star clusters with ages between 

1 and 40 Myr, and high present star-formation rates (SFR) 
of ~22 M /yr ([Adamo et aLl|2010[). The ste llar mass was 
estimated as 10 ?u Mq by |Ostlni^aLl(|200l|) whereas the 
total gas mass was derived by Bergvall et al.| pOOO) to be 

2 x 10 9 M©. An intriguing aspect of Haro 11 is that it is a 
Lya emitting galaxy and can be considered as both a lo- 
cal (albeit faint) analogue of the high redshift Lyman Break 
Galaxies (LBGs) or a rather bright Lyman Alpha Emitter 
(IHayes et al. 2007 Leitet et al. 20111. Moreover, it is also 



in the rare class of objects that are Lyman continuum emit- 
ters and its IR luminosity of ~ 2 x 10 11 L© classes it as a 
luminous IR galaxy (LIRG) ( | Adamo et al.|2010 l. 

This paper focuses on another interesting aspect of 
Haro 11 - its classification as a 'high N/O' galaxy. The 
distribution of nitrogen as a function of metallicity in star- 
forming galaxies (SFGs) has often been a cause for debate, 
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especially at intermediate metallicities (7.6^ 12+log(0/H) 
^8.3). Within this metallicity range the N/O abundance ra- 
tio indicates that nitrogen behaves neither as a primary ele- 
ment (i.e. N/O independent of O/H) nor as a secondary ele- 
ment (i.e. N/O proportional to O/H); instead a large scatter 
is observed (see e.g. Lopez-Sanchez & Esteban 2010b, their 
fig. 11). With a nitrogen excess factor of ~6 for its metallic- 
ity having been reported for Haro 11 ( Izotov fc Thuan|1999[ ), 
this galaxy belongs to a small subset of high N/O galaxies 
highlighted by Pustilnik et al.| ( |2004[ ) , several of which have 
been the focus of our integral-field spectroscopic (IFS) ob- 
servations ( James et al.||2009| |2010| |2013 1 , aimed at under- 



Table 1. General Properties of Haro 11 



standing the source of such abundance peculiarities. In the 
past, long-slit chemical abundance studies have suggested 
that a connection exists between anomalously high N/O ra- 
tios and WR stars, based on their simultaneous detection. 
However, spatially resolved abundance studies are showing 



that this is not a one-to-one relationship (e.g. Monreal-Ibero 



et al. 2012|. IFS observations of high N/O galaxies have 



distinct advantages. By mapping the electron temperatures 
(T e ) one avoids measuring integrated and therefore likely in- 
appropriately flux-weighted line ratios, enabling more accu- 



rate chemical abundances to be computed (e.g. Kobulnicky 



et al. 1999 James et al. 2010 2013 Pilyugin et al. 20121 



Secondly, the data allows one to isolate regions of enrich- 
ment, and thirdly, these regions can be related to environ- 
mental factors such as WR-populations, large scale outflows, 
mergers etc. ( Walsh fc Roy|1989 Lopez-Sanchez et al.|2007 



James et al.|2009||Perez-Montero et al.|2011||Monreal-Ibero| 



et 



al.|2012 i 



Haro 11 is known to be a complex system. Its dynam- 
ics base on the Ha emission have been studied by 



Ostlin 



et al. (2001 1, who reported a multi-component velocity field 



and a perturbed morphology. These have been interpreted as 
signatures of a merger be tween a low-mass, evolve d system 
and a gas-rich component. Bergvall & Ostlin ( 2002 1 detected 



strong emission from WR stars in the spectral region around 
He II A4686 that suggests that Haro 11 may harbour a pop- 
ulation of such stars, although to date this population has 
not been studied in detail. 

In this paper we present integral field spectroscopy ob- 
servations which afford us a new spatiokinematic '3-D' view 
of Haro 11. The aims are three-fold; firstly to map the mul- 
tiple gas components, secondly to perform a chemodynami- 
cal study, i.e. determine the range of physical properties and 
chemical compositions for each emission line velocity compo- 
nent (e.g. Esteban fc Vilchez|l 992 Jam es et al.|20 09 Hagcle 



et al.|2012 Amorin et al.|2012 1, and thirdly, to spatially re- 
late any chemical anomalies to environmental factors. The 
adopted distance for Haro 11, along with several other of its 
properties are listed in Table [I] 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of Haro 11 were obtained with the Fibre Large 



Array Multi Element Spectrograph, FLAMES ( Pasquini et al. 
[2002] ) at Kueyen, Telescope Unit 2 of the 8.2 m VLT at 
ESO's Paranal observatory, in service mode on the dates 
specified in Table [2] Observations were made with the Ar- 
gus Integral Field Unit (IFU), with a field-of-view (FoV) 
of ll'.'5x7'.'3 and a sampling of O'.'52/lens. In addition to 



Property 



Value 



Reference 



RA (J2000) 
Dec (J2000) 

z 

Distance / Mpc 
Velocity / kms -1 
Stellar mass / Mq 
Gas mass / Mq 
Cluster masses / Mq 
Present SFR / M yr -1 
Lir I Lq 



00h36'52.5" 
-33deg33'19" 
0.020558 

83.6 
6146±17 
~10 10 
2 x 10 9 
10 3 -10 6 
22±3 
1.9 x 10 11 



NED a 
NED a 
this work 
this work b 
this work 



Adamo et al. 



Bergvall et 
|Adamo et aTj ^OlO ) 
Adamo et al. (2010) 
Bergvall et^T| |200(ty 



a NASA/IPAC Extragalactic Database 
" Derived for a Hubble constant of H 
|Bernardis et al.|2008) . 



science fibres, Argus has 15 sky-dedicated fibres that simul- 
taneously observe the sky and that were arranged to sur- 
round the IFU FoV. The positioning of the IFU aperture 
for Haro 11 is shown overlaid on tfST-ACS/WFC F435W 
imaging in Fig. [l] 

Four different low-resolution (LR) grisms were utilised: 
LR1 (AA3620-4081, 24.7±0.2 kms" 1 FWHM resolution), 
LR2 (AA3964-4567, 24.9±0.1 kms" 1 ), LR3 (AA4501-5078, 
24.7±0.2 kms" 1 ) and LR6 (AA6438-7184, 22.2±0.3 kms" 1 ). 
This enabled us to cover all the important emission lines 
needed for an optical abundance analysis, as detailed in Sec- 
tion [3] The seeing (0.45-0.79"), airmass and exposure time 
for each data set are detailed in Table [2] In addition to the 
science frames, continuum and ThAr arc lamp exposures as 
well as spectrophotometric standard star observations were 
obtained (also detailed in Table [2b. 

The data were reduced using the girBLDRS pipeline 
(Blecha & Simond 2004^ in a process which included bias 
removal, localisation of fibres on the fiat-field exposures, ex- 
traction of individual fibres, wavelength calibration and re- 
binning of the ThAr lamp exposures, and the full processing 
of science frames which resulted in flat-fielded, wavelength- 
rebinned spectra. The pipeline recipes 'biasMast', 'locMast', 
'wcalMast', and 'extract' were used (cf. |Tsamis et al.|2008| >. 
The frames were then averaged using the 'imcombine' task 
of IRAlj^] which also performed the cosmic ray rejection. 
The flux calibration was performed within IRAF using the 
tasks calibrate and standard. Spectra of the spectropho- 
tometric stars quoted in Table[2]were individually extracted 
with girBLDRS and the spaxels containing the stellar emis- 
sion were summed up to form a ID spectrum. The sensitiv- 
ity function was determined using SENSFUNC and this was 
subsequently applied to the combined Haro 11 science ex- 
posures. The sky subtraction was performed by averaging 
the spectra recorded by the sky fibres and subtracting this 
spectrum from that of each spaxel in the IFU. Custom-made 

1 The girBLDRS pipeline is provided by the Geneva Observatory 
http:/ /girbldrs. sourceforge.net 

" IRAK is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy 
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A Chemodynamical Study of Haro 11 3 



5"~ 1.8 kpc 




Figure 1. HST colour composite image of Haro 11 (consisting of ACS-HRC F220W, F330W and F814W along with ACS-WFC F435W 
and F550W images) overlaid with the 11" 5 X 7'.'3 FLAMES IFU aperture (dashed line). Knots are labeled following the nomenclature 
of |Kunth et al.| l |2003| . 



Table 2. FLAMES Argus observing log 



Date Grism Wavelength Range Exp. time Avg. Airmass FWHM seeing Standard Star 







(A) 


(s) 




(arcsec) 




03/07/2009 


LR1 


3620-4081 


3x218 


1.03 


0.45 


Feige56 




LR2 


3964-4567 


4x235 


1.02 


0.55 


LTT7987 




LR3 


4501-5078 


3x218 


1.02 


0.48 


Feige56 


19/07/2009 


LR6 


6438-7184 


4x600 


1.02 


0.64 


FeigellO 


20/07/2009 


LR1 


3620-4081 


4x218 


1.16 


0.71 


HR5501 




LR2 


3964-4567 


8x235 


1.13 


0.59 


FeigellO 




LR3 


4501-5078 


6x218 


1.06 


0.67 


FeigellO 


23/07/2009 


LR1 


3620-4081 


9x218 


1.08 


0.59 


FeigellO 




LR2 


3964-4567 


8x235 


1.05 


0.60 


FeigellO 




LR3 


4501-5078 


6x218 


1.04 


0.68 


FeigellO 


24/07/2009 


LR6 


6438-7184 


4x600 


1.08 


0.79 


LTT1020 



scripts were then used to convert the row by row stacked, 
processed CCD spectra to data cubes. This resulted in one 
science cube per grating, i.e. four science cubes for the source 
galaxy. 



Observations of an object's spectrum through the 
Earth's atmosphere are subject to refraction as a function 
of wavelength, known as differential atmospheric refraction 
(DAR). The direction of DAR is along the parallactic an- 
gle at which the observation is made. Following the method 
described in James et al. (2009), each reduced data cube 
was corrected for this effect using the algorithm created by 



3 MAPPING OF LINE FLUXES AND 
KINEMATICS 

The distribution of H I Balmer line emission across Haro 11, 
which is indicative of current massive star formation activ- 
ity, was used to delineate 3 emission line regions surrounding 
the main knots of emission apparent in the HST ACS im- 
age (Fig. [y. Our analysis of Haro 11 is based on properties 
averaged over each of these three areas whose boundaries 
are displayed on the Hcv map in Figure [2] and are labeled 
as Knots A, B, and C following the nomenclature of |Kunth] 
ennp003|. 



Walsh & Roy ( 1990 1 



Knots B and C on Fig. [T] correspond to unresolved 
super-star-clusters, with Knot B containing the youngest 
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Figure 2. H/3 emission per 0''52x0'.'52 arcsec 2 spaxel. Left - Three regions were chosen across Haro 11 whose integrated spectra were 
extracted for analysis (labelled as Knots A, B and C, after the nomenclature of |Kunth et al.| j2003|). Right - The Ha line profile 
integrated over each region is shown along with a three-component Gaussian fit. North is up and east is to the left. 



cluster population of mass > 8 x 10 6 Mq (Adamo et al. 
2010). Knot C is the brightest UV source in Haro 11 proba- 
bly due to ongoing massive star formation activity (Hayes et 
al. 2007). Knot A on the other hand is resolved to a sparse 
population of clusters. Both Knots B and C are ringed by 
dust lanes and Knot B appears almost fully surrounded by a 
dusty spiral arm inclined to the plane of the sky. Not only the 
surface brightness, but also the composite structure of the 
lines varies considerably across Knots A, B, and C (Fig. [2). 

The high spectral resolution and signal-to-noise (S/N) 
ratio of the data allow the identification of multiple velocity 
components in the majority of emission lines seen in the 300 
spectra of Haro 11 across the IFU aperture. Following the 
automated line-fitting procedures outlined in | James et al.| 
(2013 and references therein) and likelihood ratio methods 



described by Westmoquette et al. (2011 their Appendix A), 



we were able to find the optimum number of Gaussians re- 
quired to fit each observed profile. S/N ratio maps (used 
for making S /N cuts) were first made for each emission line, 
taking the ratio of the integrated intensity of the line to that 
of the error-array produced by the FLAMES pipeline^] on a 
spaxel-by-spaxel basis. Single or multiple Gaussian profiles 
were then fitted to the emission lines, restricting the mini- 
mum FWHM to be the instrumental width (see Section |2J). 



Suitable wavelength limits were defined for each line and 
continuum level fit. 

Most of the high S/N emission line profiles (i.e. all 
Balmer lines, [O n] AA3727, 3729, [O III] A495^][N n] A6584, 
[S n] AA 6716, 6731) were optimally fitted with a narrow 
Gaussian (FWHM<200 kms -1 , component CI hereafter), 
an underlying broad Gaussian component (FWHM 200- 
300 km s~ 1 , C2) , and a third narrow Gaussian (C3 hereafter) 
which can appear either red- or blue-shifted with respect to 
CI. With respect to the properties of each Gaussian com- 
ponent (i.e. FWHM and/or velocity), we consider them to 
arise from gas of different physical conditions. Towards the 
outer galactic regions the lines are optimally fitted with only 
the CI and C2 components. A few lines (e.g. He I A4026) 
are optimally fitted with a single, narrow component (CI). 
Figure [4] shows the complex profile of the Ha line, along 
with maps of its constituent velocity components. 

Tables lAll IA3I list measured FWHMs and observed and 
de-reddened fluxes for the fitted Gaussian components of the 
detected emission lines in the three regions defined in Fig. [2] 
The fluxes are for summed spectra over each SF knot and are 
quoted relative to the flux of the corresponding H/3 compo- 
nent. In addition, table|3]lists the observed and de-reddened 
integrated line fluxes (i.e. C1+C2+C3) from summed spec- 
tra over the entire galaxy. All fluxes were corrected for red- 



3 The uncertainties output by the pipeline include uncertainties 
from the fiat-fielding and the geometric correction of the spectra 
and not just photon-noise. 



4 [O III] A5007 is not within the wavelength range of this data on 
account of its redshift. 
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Figure 3. Map of E(B-V) for the integrated flux of the emission lines. Ha contours are overlaid in white. North is up and east is to the 
left. 
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14 


H8+He I A3888 
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81±1 


40 


22.28±1 


83 


He I A4025 


1 


59±0 
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1.70±0 


14 


[S II] A4068 


2 


05±0 


11 


2.19±0 


15 


HS 


23 


21±0 


48 


24.51±1 


22 


H7 


43 


63±0 


92 


45.32±2 


18 


[O III] A4363 
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31±0 


18 


2.40±0 


21 


He 1 A4471 


3 


55±0 


11 


3.65±0 


19 


[Fe 11] A4658 
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76±0 
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17 


H,3 


100 
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80 
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99 


41.22±1 


40 


He 1 A6678 


2 
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07 
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25 


28±0 
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13 


± 0.02 




F(H/3) X10 14 ergs -1 cm -2 




37.8078 ± 0.77 





Table 3. Haro 11 fluxes and de-reddenened line intensities (rel- 
ative H/3 = 100) for the integrated spectrum of Haro 11. Line 
fluxes were extinction-corrected using the c(H/9) value shown at 
the bottom of the table, calculated from the relative Ha, H/3 and 
H7 fluxes. The line intensities listed here were used for summed- 
spectra T c and N c diagnostics and ionic abundance calculations 
listed in Table g] 



dening using the galactic reddening law of Howarth ( 1983 1 
using c(H/3) values derived from the Ha/H/3 and H7/H/3 
line ratios of their corresponding components, weighted in a 
3:1 ratio, respectively, in conjunction with theoretical Case 



B ratios from Storey & Hummer (1995). A Milky Way red- 
dening of E(B — V) = 0.011 mag, in the direction to Haro 11 



is indicated by the extinction maps of Schlegel et al. ( 1998 1, 
corresponding to c(H/3)=0.016. Total c(H/3) values applica- 
ble to the galaxy and its individual emission line components 
in each main SF knot are listed in their respective tables. 

An extinction map was also derived using the Ha/H/3 
and H7/H/3 emission line ratio maps and is shown in Fig. [3] 



This was computed using the integrated line fluxes (i.e. 
C1+C2+C3). The distribution in E( B-V) is in agreemen t 
with that found in the Lya study by Ha yes et al.| ( |2007[ ), 
where the extinction appears to peak at the locations of 
Knot B and C whereas Knot A shows little or no extinc- 
tion. As pointed out by |Hayes et al.| ( |2007[ ), it is interesting 
that the source of Lya emission within Haro 11 is located in 
Knot C, whereas one would expect a low level of extinction 
to be needed to allow the escape of Lya photons. This indi- 
cates that dust is not the major regulatory factor governing 
the Lya morphology in Haro 11, which was found by [Hayes] 
|et al. I ( |2007[ ) to be driven more by the H I and kinematical 
structure. This is similar to what is found in many other 



Lya emitting galaxies (e.g. Kunth et al. 1998 Mas-Hesse 



|et al. 2003 ; H ayes et al.||2010| and references therein) 



3.1 Ha maps 

In the light of Ha (Fig. [4^l) Haro 11 consists of three 
main regions: a central body located at Knot B (a = 
00 h 36 m 52.41 ± 0.1 s , 6 = -33° 33' 16.82 ± 0.1", J2000), and 
two separated peaks, Knot C (~3'.'l, ~1.2 kpc to the east) 
and Knot A (~3'/2, ~1.3 kpc to the south). Knots A and B 
are clearly evident in all 3 line components, with Knot B ap- 
pearing most pronounced in C2. Knot C appears as a lower 
surface brightness feature in CI and C2, projected against 
a fainter background connecting it to Knot B. The inter- 
knot background appears fainter in C3. In integrated light 
(C1+C2+C3), the Ha flux of ~ 1.06 x 10~ 12 ergs cm -2 s~\ 
in reasonable agreement with 



et al. (19931 



Ostlin et al. 



(1999| and 



Vader 



The Ha radial velocity maps of Haro 11 are shown in 
Fig. [4j> (relative to the systemic velocity). The CI veloc- 
ity component shows a velocity gradient orientated in the 
NE-SW direction, ranging from ~ —50 to +100 kms -1 , and 
an increase up to ~200 kms -1 along the western edge. A 
tentative rotation axis can be defined based on component 
CI which is at an angle of 90 deg to the dusty spiral lane 
surrounding Knot B in Fig. 1. In contrast to this, the veloc- 
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Figure 4. Maps of Haro 11 in the Ho velocity components C1-C3. The first row shows an example three-component Ha line pro- 
file. The highlighted component indicates to which component the maps in each column correspond; a) logarithmic flux (in units of 
10 -16 ergs s~ 1 cm~ 2 arcsec~ 2 ); b) radial velocity (relative to the heliocentric systemic velocity of 6146 kms" 1 ); c) FWHM corrected 
for the instrumental PSF. See text for details. North is up and east is to the left. 



ity gradient of C2 ranges from ~ —20 to +100 kms~ 1 and 
extends in the NW-SE direction across Knot B, i.e. orthog- 
onal to that of CI. This peculiar v elocity structure is also 
highlighted by Ostlin et al. (20011, who suggest that this 



central region contains 'a counter-rotating disc or high ve- 
locity blobs', which we are tracing with component C2 and 
is centred on Knot B. This counter-rotating disc may also 
be evident as dusty filaments to the east of Knot B cut- 
ting across the spiral that is ringing it (Fig 1). Overall, and 
in accordance with the conclusions of Ostlin et al. (2001), 



the combined velocity distributions of CI and C2 suggest 
that the center of Haro 11 is not dynamically relaxed, while 
the outer velocity field shows a rather slow rotation. Com- 
ponent C3 shows a more patchy velocity distribution but 
there is evidence for a wide-angled bipolar outflow oriented 
approximately along the north-south direction centred on 
Knot B. The overall velocity distribution of C3 could also 
be interpreted as the lower S/N ratio signature of geomet- 
rically thin shells created by stellar outflows and supernova 
remnants projected against the background emission. 
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Figure 6. De-reddened emission line maps of Haro 11 in 
[O II] A3726/A3729, with integrated light Ha emission line con- 
tours overlaid in dashed line. North is up and east is to the left. 



in Fig. [5}i, in the CI and C2 velocity components. Unlike 
the CI components of the [N n] and [O in] emission lines, 
[S n] A6716 does not show clear peaks in surface brightness, 
but instead remains relatively constant throughout the main 
body of Haro 11. In comparison, the C2 velocity component 
of [S n] A6716 peaks strongly in Knot B with a morphology 
akin to Ha C2 emission (Fig. [4^, - central panel). Figure [5]b 
shows Haro 11 in the light of [Ne ill] A3868, which was only 
present in a narrow (CI) component, that has three clear 
peaks at the three main knots of star-formation. 

Figure [5] shows Haro 11 in the light of the electron den- 
sity (JV e ) sensitive doublet ratio, [O n] AA3726/3729. A small 
amount of structure is visible; the ratio peaking at the lo- 
cation of Knot B, which suggests a region of increased N e . 
The broad component of the [O il] AA3726, 3729 doublet 
cannot be reliably decomposed as their relatively small rest 
wavelength separation is blurred out by its large width. The 
[S n] AA6716, 6731 doublet on the other hand is unusable as 
at the redshift of Haro 11 the latter line falls in a spectral 
region of strong telluric oxygen absorption. 



FUSE observations of Haro 11 by Grimes et al. (20071 



report similar findings to the Ha velocity maps presented 
here, albeit in one-dimension and relative to a slightly higher 
systemic velocity of 6180 kms -1 . UV absorption lines show 
strong absorption features (with FWHMs of ~300 kms -1 ) 
blue-shifted by 100 kms -1 relative to the host galaxy sys- 
temic velocity (rather equivalent to our C2 component) and 
an additional, weaker high-velocity outflow +200-280 kms -1 
away from the galaxy (our CI component). Interestingly, 
when comparing the equivalent widths, the higher ionization 
UV lines are stronger in the high-velocity outflow than the 
low-velocity absorption features, which Grimes et al. (2007) 
suggest could be due to shocks between the outflow and the 
ambient material. 

The maps of Ha FWHM (Fig. EJ;) show minimal struc- 
ture in each of the three main velocity components. FWHMs 
range from ~80-180 kms -1 for CI, ~180-450 kms" 1 for 
C2 and <90 kms -1 for C3. |6stlin et ak] |200l| observe a 
maximum FWHM of ~270 kms -1 for the emission lines 
within Haro 11, which agrees well with the average FWHM 
of 283±5 kms _1 seen for C2 throughout the main body of 
the galaxy. The largest FWHM value for C2 is observed di- 
rectly east of Knot A coinciding with a super star cluster 
positioned at the end of a bright arc of emission extend- 
ing southwards from Knot B (Fig. 1). This region spatially 
aligns with a region of blueshifted C3 emission and could 
thus be a further signature of outflow activity originating in 
Knot B. 



3.2 Forbidden Lines 

Emission line maps of Haro 11 in the light of the [O ill], 
[N n], [S n], and [Ne ill] lines are shown in Fig. [5] The 
morphology of [O ill] A4959 is similar to that of Ha, with 
the peak in C2 emission being more extended than that of 
CI. Despite the low S/N of the [O in] A4363 line, we were 
able to map the CI velocity component throughout Haro 11, 
with peaks that correlate with those in Ha along with the 
C2 component within Knot B. 

We show Haro 11 in the light of [S n] A6716 emission 



4 DIAGNOSTICS: ELECTRON 

TEMPERATURE AND DENSITY 

The dereddened [O III] (A4959)/A4363 and [O n] 
A3726/A3729 intensity ratios were used to determine 
electron temperatures and electron densities throughout 
Haro 11, respectively. The T c and N c values were computed 
by inputting the [O ill], and [O n] intensity ratios into 
IRAF's ZONES task in the NEBULA package. This task derives 
both the temperature and the density by making simulta- 
neous use of temperature- and density-sensitive line ratios. 
Atomic transition probabilities for 2+ and + were taken 
from 1 
from 



Wiese et al. ( 1996[) w hilst collision strengths were taken 
LennonfeBurke|(l994| ) and |McLaughlin fc Bell| ( p93l ), 



respectively. 

T e maps computed for line components CI and C2 are 
shown in Fig. [7] and the N c map for CI is shown in Fig. [8] 
The per spaxel S/N ratio of the T e maps is in the range 10- 
20. The error- weighted mean values for each of the Knots 
(i.e. map-based averages across apertures defined in Fig. [2| 
are listed in Table [4] along with uncertainties correspond- 
ing to the standard-deviation across each aperture. Fig. [7[i 
shows evidence for a small amount of T c variation through- 
out Haro 11. For component CI Knots A and B show tem- 
peratures of ~12-13,000 K whilst Knot C shows ~16,000 K. 
The broad line component in Knot B has a lower tempera- 
ture of ~9,000 K. 

Electron density variations are also apparent and denser 
areas are correlated with Knots A and B which show den- 
sities of ~400 cm' 3 whereas Knot C is less dense with 
~150 cm -3 . The low N e seen in Knot C may help ex- 



plain the escape of Lya photons within this region (Hayes 
|et al.|p007| . In Knot C, |Hayes etH] p007) found that the 
Lya/Ha ratio agrees with the theoretical case B recombina- 
tion value, Lya/Ha =8.7 (Osterbrock & Ferland 2006), indi- 
cating that the Lya photons from this region are scattered 
just a few times before escaping. Such direct escape would 
be possible from a low-density, dust-poor nebula bordering 
on the optically thin Case A, however, as we discussed in 
the introduction to Section 3, the dust content is probably 
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Figure 5. Dereddened intensity maps of Haro 11 in [O III] A4959, [N n] A6584, [O III] A4363 and [Ne III] A3868 for the separate velocity 
components CI, C2 and C3. [O III] A4363 is shown after being re-binned by 1.5x1.5 spaxels and re-mapped onto the original grid (see the 
text for details). Only spaxels with S/N^10 are shown for [O III] A4959 and [N II] A6584, ^5 for [Ne III] A3868 and >1 for [O III] A4363. 
North is up and east is to the left. 



not the determinant factor in this case; the lower density 
may be the cause. As suggested by the models of |Neufe"Id| 
(19911, Lycv photons will suffer less attenuation than radi- 



ation that is not resonantly scattered when they exist in a 
multi-phase medium, i.e. dusty gas clouds embedded within 
a low-density inter-cloud medium. Our N c values based on 
the CI [O ii] A3726/A3729 ratio is not in agreement with that 
of Guseva et al. (20121; this may be because their [S n] di- 



agnostic is affected by telluric absorption of the A6731 com- 
ponent. 



5 CHEMICAL ABUNDANCES 

Ionic abundance maps relative to H + were created for 
the N + , 0+, 2+ , Ne+ and S + ions, using the AA6584, 
3727+3729, 4959, 3868 and 6717 lines respectively. Exam- 
ples of flux maps used in the derivations of these ionic abun- 
dance maps can be seen in Figure [5] Abundances were cal- 
culated using the ZONES task in IRAF, using the respec- 
tive T e and N e maps described above, with each FLAMES 
spaxel treated as a distinct 'nebular zone' with its own set 
of physical conditions. The IRAF abundance results are cor- 
roborated by comparing with abundance maps computed via 
the programme equib (originally written by I. D. Howarth 
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Figure 7. Electron temperature map of Haro 11 for velocity components CI and C2. Overlaid circles represent the size of the uncertainty 
in T c within each spaxel. North is up and east is to the left. 




tionship between the S 2+ and S + ionic fractions from the 
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Figure 8. Electron density map of Haro 11 for velocity compo- 
nents CI and C2. Overlaid circles represent the size of the uncer- 
tainty in T e within each spaxel. North is up and east is to the 
left. 



and S. Adams), which uses different atomic data tables; we 
found that typically the differences in the abundance ratios 
are within 5%. 

Ionic nitrogen, neon and sulphur abundances were con- 
verted into N/H, Ne/H and S/H abundances using ioni- 
sation correction factors (ICFs) from Kingsbu rgh &: Bar- 
|low| (jT994) . The O/H abundance was obtained by adding 
the O^/H 4 " and + /H + abundance maps (i.e. assuming 
that the 3+ /H+ ratio is negligible). Since the [S ill] A6312 
line was not within the wavelength range of our data, the 
S 2+ /H + abundance was estimated using the empirical rela- 



corrected equation A38 of Kingsburgh & Barlow (19941, 
namely S 2+ /S+ = 4.677 x (CH + /0 + ) u "' M . The S/H results 
should therefore be considered with caution, but they are not 
at variance with those obtained using the ICF prescriptions 
of Izotov et al. (2006). O/H, N/H and log(N/0) abundance 
maps for line component CI are shown in Fig. [9] and arc 
described below. Abundances for the three regions were de- 
rived from taking error-weighted averages across the regions 
in the abundance maps shown in Fig. [9] These map-based 
spatial averages are listed in Tableland correspond to the 
'average values' described the following section. We were 
unable to derive ICFs for the C2 line component without 
maps of [O n]-C2, however, averages over the ionic abun- 
dance maps are listed in Table [4] 

5.1 Maps 

The elemental oxygen abundance map shown in the top 
panel of Fig. [5] displays a single peak whose location cor- 
relates with Knot B (the brightest of the knots in terms of 
Ha flux). Both Knots A and C reveal a decreased oxygen 
abundance at the location of their peak in Ha flux. Slight 
abundance variations are seen across Haro 11, with a de- 
creasing oxygen abundance as one moves away from Knot 
B, with a average metallicity of 12+log(0/H)= 8.25±0.03, 
decreasing to 8.09±0.03 and 7.80±0.04 for Knots A and C, 
respectively. There is, therefore, a distinct inhomogeneity in 
oxygen abundance across Haro 11, with Knot B having a 
higher oxygen abundance by upto ~0.4 dex. Whilst the re- 
sult for Knot B is in relatively good agreement with |Guseva| 
et al.| (|2012[> who derived 12+log(Q/H)= 8.33, our abun- 



dance for Knot C is somewhat lower than their 8.10 due to 
our higher temperature (they derive ~1 1,500 K for Knot C). 
We are, however, in agreement with the va lue 12+log(0/H) 
= 7.9 obtained bylBergvall & 6stlin||2002| in a 4" x 4" aper- 



ture over the brightest region of Haro 11 in the visual region. 
Whilst we cannot comment on the C2 elemental abundance 
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Figure 9. Elemental nitrogen, oxygen and N/O abundance maps 
for velocity component CI, with integrated light Ha emission line 
contours overlaid in dashed line. North is up and east is to the 
left. 



of oxygen, Table[4]shows that the 2+ C2 component abun- 
dance in Knot B is approximately twice that of its narrow 
CI component. 

The middle panel of Fig. [9] shows the N/H abundance 
across Haro 11. The peak in abundance is located within 
Knot B, with an average N/H ratio of 9.38 ± 0.02 x 10" 6 , 
which decreases to ~ 4.88 x 10~ 6 in Knots A and C. In 



contrast to oxygen, there is no distinct decrease in nitrogen 
abundance in Knot C and as a result this knot shows a 
peak in N/O, which can be seen in the bottom panel of 
Fig. [9] Whilst Knots A and B display relatively 'normal' 
log(N/0) levels for their respective metallicities (~-1.4 and 
-1.3, |L6pe z-Sanchc z &, Esteban| (see |2010b their figure 11)), 
Knot C has an average log(N/0) of -1.12, ^0.5 dex higher 
than the expected value for 12+log(O/H)=7.80. In fact, the 
map of log(N/0) shows that this peak actually reaches ~- 
1.0 to —0.8 at the boundary between Knots B and C (as 
can be seen by the Ha contours). Therefore, the factor of 
~3 excess reported by |Pustilnik et al.| ( 2004 1 is seen by us 
within Knot C only. We investigate the cause behind this 
anomalously high ratio in Section [6] 

Again, we cannot draw firm conclusions on the N/H or 
N/O abundance within the broad velocity component gas. 
However, as a rough estimate, if we apply the average ratio 
of + /0 2+ = 1.57 seen in the CI component of Knot B to the 
2+ /H+ value in C2, we would obtain 0+/H+ = 2.43 x 10" 4 
and a log(A r+ /0 + ) ratio of —1.07. This would provide some 
evidence for the existence of a broad-line region within Knot 
B that has a higher O/H abundance than the narrow-line 
component, but a similar N/O ratio. 

Guseva et ah] poll} derive log(N/0) values of -0.92 



and —0.79 for Knots B and C, respectively; whilst the value 
for Knot C is in agreement with the maps derived here, 
the value for Knot B is somewhat higher. We interpret this 
as being a result of their using ratios from integrated line 
profiles. 

The Ne/H distribution in Haro 11 (not shown) follows 
that of nitrogen in that a peak in abundance is seen in Knot 
B with Ne/H = l.OxlO" 4 , whilst Knots A and C have an 
abundance ratio of ~ 0.4 x 10 -4 . The average Ne/O ratios 
are all ~0.4 dex higher than the expected range of log(Ne/0) 
~ -0.9 to -0.4 ( |L6pez-Sanchez fc Esteban|2010b"| . Deriva- 
tions of Ne/H were also made using the neon ICF of |Izotovj 



et al. (20061, with only the Ne/O ratio for Knot B being 



-^0.1 dex higher than the expected value. The S/H distribu- 
tion across Haro 11 (also not shown) is relatively constant 
for the CI gas component, with each region displaying aver- 
ages within the range of ~(3-6) x 10~ 7 . The S/O ratio for 
Knot C is in agreement with the average reported range of 
log(S/0> 1.9 to -1.5 ( jlzotov et aLl|2006[ > for blue com- 
pact galaxies, whilst Knots A and B are ~0.2-0.4 dex lower 
than the expected values. 



5.2 Comparison of spatially resolved and global 
spectra abundances 

Also listed in Tableware the T c , N c , ionic and elemental 
abundances derived from the galaxy integrated spectra (us- 
ing the line fluxes listed in Table [3j|. Following the method- 
ology of James et al. (20131, this set of results allows us to 



assess if integrated galaxy spectra can reliably represent the 
physical properties of the ionized ISM. In a comparison be- 
tween spatially-averaged and integrated-spectra results for 



UM 448, James et al. (2013) found that whilst relative abun- 



dance ratios of heavy elements can be reasonably well ob- 
tained from integrated spectra and/or long-slit observations, 
the ionic and elemental abundances relative to hydrogen are 
not well reproduced. |James et al.| |2013) suggested that this 
may be because the former properties do not have a very 



A Chemodynamical Study of Haro 11 11 



strong sensitivity to T c biases and/or, in contrast to inte- 
gral field spectroscopy, integrated properties derived from 
flux-averaged methods cannot resolve abundance and tem- 
perature variations and sites of enrichment are being aver- 
aged out. 

The present results suggest that whilst the N/O and 
Ne/O values are in agreement with the spatially-averaged 
values, the ionic and elemental abundance ratios relative to 
hydrogen are not. The T e from the integrated spectrum is in 
excellent agreement with the A4959-flux-weighted mean T e 
across the three regions (f0,850 K), suggesting that the in- 
tegrated [O ill] ratio may have been more heavily weighted 
by areas of stronger [O ill] A4959 emission. Indeed, the aver- 
age T c for Knot B (which has the strongest [O in] emission) 
from the CI and C2 maps, would be ~10,550 K - in agree- 
ment with the integrated-spectrum T c . These discrepancies 
suggest that values derived from long-slit observations may 
be prone to luminosity-weighted fluxes and resultant inaccu- 
rate ratios, compounded by the fact that large apertures can 
include a mixture of gas with different ionization conditions 
and metal content. Such effects can lead to significant biases 
in analyses based on integrated spectra of galaxies both in 
the local and the high-redshift universe. These biases may 
be of comparable magnitude to the effects of temperature 
or density variations (e.g. Peimbert & Costero 1969; Rubin 



1989 Viegas &: Clegg||l994 | often invoked in nebular abun- 
dance analyses. 



5.3 Helium Abundances 

The He + /H + abundances for Knots A and B are listed in 
Table [6] The He I lines in Knot C are abnormally weak rel- 
ative to H/3 probably due to the effects of underlying stellar 
absorption and hence a reliable He abundance cannot be 
derived for this area. Summed spectra over each knot and 
velocity component (i.e. C1+C2-I-C3) were used rather than 
map-based averages because of the low per spaxel S/N ra- 
tio of the lines. Fluxes were de-reddened using the c(H/3) 
value for its corresponding knot. Helium abundances were 
calculated using the error-weighted line fluxes listed in Ta- 
ble 5 and adopting the map-based average T c and N e for 
each knot (listed in Table |4j. The Case B He I emissivities 
of Smits (1996) were adopted, correcting for the effects of 



collisio nal excitation using the formulae in [Benjamin et al"] 



I 1999). The presence of neutral helium within the H~^zone 



has b een accounted for using the ICF of Pei mbert et al.| 
( 1992). No substantial amounts of He 2+ are expected to be 



present. 



6 WR STARS AND N-ENRICHMENT 

As discussed previously, a connection is thought to exist be- 
tween galaxies with high N/O ratios and the presence of 
Wolf-Rayet (WR) stars (e.g., |Brinchmann et aTl|2008[ and 
references therein). Two examples are NGC 5253 (Lopez- 



Sanchez et al. 2007| and references therein) and Mrk 996 
( ] James et al.|2009[ ), where spatially-resolved studies showed 
regions of enhanced N/O that correlate with a strong WR 
population. However, the relationship may not be clear-cut. 
Cases also exist where galaxies display N-overabundance 
over areas too large for WR-stars to pollute, suggesting that 



other mechanisms may be responsible ( Perez- Montero et al. 
2011||James et al.|2013"l ) 



The presence of WR features within the spectra of 
Haro 11 was reported by Bergvall & Ostlin (2002). However, 



to-date, no attempt has been made to map or quantify the 
size of this population. The FLAMES-IFU spectra also show 
a broad WR stellar feature around 4690 A, attributable to 
a mixture of late-type WN (WNL) stars, as seen in Fig. 10 
A map of the WR feature is shown in Fig. 11 (made by inte- 



grating the flux across the feature) , where the WR emission 
is seen to extend throughout Knots B and A, with only a 
minimal detection in Knot C. 

When searching for WR features one must consider the 
spatial width of the extraction aperture, which can some- 
times be too large and thereby dilute weak WR features 
with continuu m flux (|Kehrig et al.||2008| |L6pez-Sanchez fc| 
|Esteban|20"08||James et al.|2010||Perez-Montero et al.|2011[ ). 
We therefore also examined the spectra summed over each 
individual knot (Fig.[5|, which are shown in Fig. |10| In order 
to determine the number of WR stars from the strength of 
the 4640-4690 A feature for the summed spectra over Knots 
A, B and C, we followed the procedure described in James et 
al. (2012). This ma kes use of an LMC WR (WN 5-6) spec- 
tral template from Crowther & Hadfield (2006) to fit the 



observed WR broad feature using Monte Carlo to estimate 
the errors. 

It is thus estimated that 900 ± 400, 1500 ± 300 and 
300 ± 400 WN stars are present in Knots A, B and C re- 
spectively. A total population of 5200 ± 1500 WN stars is 
estimated applying the same procedure for the integrated 
spectrum of the entire galaxy. The difference between this 
and the total number of WR stars from the regional spectra 
suggests that there may be WR stars located outside knots 
of star- format ion, similar to the situation in the nearby BCG 



NGC 5253 (e.g. |Monreal-Ibero et al.||2010| and references 
therein). 

Due to the short WR phase of a star's evolution, an- 
other way to assess if they are present is to estimate the 
current age of the ionizing stellar population. This can be 
achieved by assessing the equivalent-width (EW) of H/3 at 
a given metallicity. A map of EW(R/3) is shown in Fig. 12 



Two regions of increased EW are present, located at Knots 
B and A. Knot C, on the other hand, shows a distinct de- 
crease in EW. Following the method outlined in jjames et al.| 
( |2010| , we can use this map, in conjunction with the metal- 
licity map described in Section |5.1| to estimate the age of 
the most recent star- forming episodes throughout Haro 11, 
by comparing the regional observed average EW(H/3) values 
with those predicted by the spectral synthesis code STAR- 
BURST99 ( |Leitherer et al.|1999[). We utilize the same model 
parameters described in James et al. (20101 and obtain av- 



erage stellar ages of 4.9±0.4, 4.3±0.5 and 7.4±0.3 Myr for 
Knots A, B and C, respectively. These ages are in relatively 



good agreement with Adamo et al. (20101 who derive ages 



of 3.9 and 9.5 Myr for Knots B and C, respectively from 
modeled spectral energy distributions of the star clusters. 



According to Leitherer et al. (19991, typically an instanta- 
neous starburst shows WR features at ages of ~3-6 Myr 
for metallicities of Z = 0.004-0.008, similar to those seen in 
Haro 11. Thus the stellar population ages of Knots A and B 
are consistent with hosting WR stars, whilst Knot C appears 
too old to do so. 
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Table 4. Ionic and elemental abundances for Haro 11, with the corresponding T e and N e used in their derivation. Columns 1-4 refer 
to abundances derived from maps, weight- averaged over regions defined in Fig. [2] All ionisation correction factors (ICFs) are taken from 
|Kingsburgh fc Barlow| \1994\ . 



Due to their high mass, the number ratio of WR stars to 
O stars within a galaxy can be used to provide constraints 
on the high-mass end of the initial mass function (IMF). 



Table 5. Summary of Knot Properties 



Following the methodology outlined in James et al. (20091, 
an estimation of the number of O-type stars within each 
knot was made using L(H/3), A'wr, and the number of 07V- 



equivalent stars (e.g. Kunth fc Sargent||1981 Osterbrock & 



Cohen||1982 1, after applying the appropriate correction fac- 
tor for the proportion of other O star subtypes, rjo(t), for 
the metallicity of each knot. The emission rate of ionising 
photons, Qg 6s , are 1.01, 3.20 and 0.34 xlO 53 photons s" 1 
for Knots A, B, and C, respectively, yielding average O-star 
populations of size ~45,900, ~101,700 and ~25,800. As a re- 
sult, we find AW/AW+o ratios of 0.020, 0.014 and 0.012 
for Knots A-C. These results are in agreement with the 



predictions of Schaerer & Vacca (1998) for Knots A and C, 



but for an IMF with a = 3.0 (see 



Lopez-Sanchez & Esteban 



2010a their figure 13), i.e. slightly steeper than the Salpeter- 
like IMF with a =2.35. However, the ratio for Knot B calls 
for an even steeper IMF. One must bear in mind that here 

as [Lop ez-Sanc hez fc| 



we are assuming A'wr 



Ny 



and, 



Esteban (2010a) point out, the inclusion of WCE stars may 
increase the iVwH/^WR+o ratio sufficiently to reconcile our 
results with the standard Salpeter-IMF. 

By comparing the map of N/O (Fig. [9]) and WR emis- 
sion (Fig. Ill we are thus faced with two intriguing situa- 
tions. Firstly, for Knot C, we have a large area of gas with a 
relatively high N/O ratio and a small WR population, and 
secondly, for Knots B and A, we have the reverse situation; 
regions with relatively low N/O ratio but large WR popu- 
lations. In the following we will discuss scenarios that may 
explain these circumstances (for guidance, we summarise the 
main properties of each knot in Table [5J. 

Knot C: The negligible amount of WR emission within 
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this region, combined with the large projected area of high 
N/O suggests that N-rich winds may not be the cause of the 
observed Alog(N/O)-~+0.5 dex excess. Whilst, this enrich- 
ment may be related to the stellar wind ejecta coming from 
WR stars in neighbouring Knot B, the ~1 kpc separation 
distance between Knot B's peak in WR emission and the 
Knot C region of enhanced N/O suggests this is unlikely. 
Instead, the excess may be due to other global processes 
occurring in Haro 11. For example, the galaxy's unrelaxed 
dynamics suggests that it is a merger of two bodies. Indeed, 
the differing age, A^, T c , and metallicity of Knot C com- 
pared to the rest of Haro 1 1 (see Table [5| suggest that it 
has undergone a separate evolutionary path. 

The decreased oxygen abundance in Knot C indicates 
that accretion of metal-poor gas might be implicated; some 
studies suggest that interacting galaxies fall ~0.2 dex be- 
low the mass-metallicity relation of normal galaxies due 
to tidally induced large-scale metal-poor gas inflow to the 
galaxies central regions (e.g. |Kewley et al.||2006| |Michel-| 
Dansac et al.|200 8 Pccpl es~et al.|2009 l. However, the stellar 
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mass of Knot C (1.36 x 10 7 M |Adamo et al.|2010[ ) and its 
map-derived average metallicity (12+log(O/H)=7.8±0.13), 
actually place it in direct agreement with the M-Z relation- 
ship at low masses (see |Zahid et al.|2012 their figure 6) 



alongside low-mass dwarf irregulars from ( Lee et al. 2006 1 



and other blue compact dwarf galaxies from (Zhao et al 



20101 



Conversely, there may have been an outflow of oxygen- 
enriched gas due to supernovae (see e.g. |van Zee fe Haynes| 
20061, rather than an inflow of metal-poor gas. The effec- 



tive yield of oxygen and nitrogen may not be in-sync since 
the former is formed in high-mass stars, and the latter in 
intermediate- mass stars, and it is the more recently enriched 
gas that is ejected via supernova winds (e.g. the case of 
NGC 1569, |Martinet al.|2002[ ). Indeed, being 2-3 Myr older 
than the rest of the galaxy, Knot C may have experienced 
more SNe explosions than the younger Knots, A and B. This 
scenario might be more plausible with the aforementioned 
top-heavy IMF that is required to explain the number ratio 
of WR/O stars. 

Alternatively, the age of Knot C (7.4 Myr) implies that 
its WR phase has recently been completed, whereas in Knots 
B and A (both with ages of ~4-5 Myr) the WR phase is 
still ongoing. This suggests that Knot C has had sufficient 
time to maximise the injection of N into the local ISM dur- 
ing its WR phase, which may still be relatively undiluted 
by subsequent mass loss of lower mass stars. This case is 
supported by Knot C having a lower current star-formation 
rate (~0.09 Mq yr _1 , as derived from the Ha flux listed 
in Tableland the Kennicutt 1998 relationship), in com- 
parison to Knot B (0.86 M yr _i ) and Knot A (0.26 Mq 
yr -1 ). The situation witnessed here, where the mixed ejecta 
have cooled but the WR-phase has finished confirms that 
the relationship between N-enrichment and the existence of 
WR-features is far from being one-to-one. In fact, Haro 11 is 
highly analogous to the nearby 'high N/O' BCD NGC 5253, 
where N-enrichment is not always directly associated with 
the p resence of WR stars (jWestmoquette et al.||2012, in| 



press Monreal-Ibero et al.|20l6" l. As more spatially-resolved 
analyses come to light on 'high N/O objects', it is becom- 
ing apparent that the simultaneous detection of both WR- 
features and enhanced N/O is most probably a function of 
starburst age and the properties of the surrounding medium. 
As proposed by Westmoquette et al. ( 2012, in press[ ), in the 
two known cases where an increase in N/O is coupled with 
a WR-popu lation, i.e. 'Knot 1' in NGC 5253 ( |Monreal-Ibero 
et al.|2012 and references therein) and the central region of 
Mrk 996 (James et al. 2009), we are only able to observe 



the enriched gas because the high density and pressure of 
the ISM surrounding these areas may have acted to impede 
the WR winds, allowing their metal-enriched ejecta to mix 
with the cooler phases and subsequently become detectable 
in the warm (~ 10 4 K) gas. 

Knots A+B: The existence of a strong WR population 
within these knots, but at normal levels of N/O suggest that 
we are observing the galaxy at an interesting point in its 
evolutionary path. This situation has also been observed in 
NGC 5253 (see results for 'Knot # 2' of |Monreal-Ibero et al.| 
2012), where WR-features were seen in a region of clusters 



~l-5 Myr in age, with normal N/O levels for its metallicity. 
Here the authors suggest that we are observing a region 
whose stars are in the process of expelling their processed 
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Figure 10. Sections of continuum-subtracted FLAMES IFU 
spectra, showing the blue 'WR bump' region within summed 
spectra over the entire galaxy and separate aperture regions 1, 
2 and 3 (see Fig. B. We detect WR emission features in the 
spectra summed over the entire galaxy and also each of the in- 
dividual aperture region spectra. Overlaid are the estimates of 
5200±1500, 900±400, 1500±300 and 300±400 WN stars, for the 
entire galaxy and Knots A, B and C, respectively, using the tem- 
plates of |Crowther fc Hadfield| l |2006| . Vertical dashed lines in- 
dicate the location of [Fe in] 4658 A and a weak, narrow, He II 
emission line at 4686 A. The dashed horizontal lines in the plots 
show the zero flux levels. 



nitrogen (which begins < 5 Myr) and have not had time 
enough for this material to cool down and mix with the 
ISM. Under the assumptions of instantaneous mixing and 
that the velocity of this gas traces the velocity under which 
the contamination of extra nitrogen propagates through the 
ISM, they estimated that the cooling down and mixing of 
expelled material can last between ~2-8 Myr over an area 
40-50 pc in diameter. For comparison, the areas of Knots A 
and B are ~0.4 and ~1 kpc 2 , respectively. If this injection of 
extra nitrogen takes place as early as ~2.5 Myr (according 
to the models of |Molla fe Terlevich||2012[ ), the N-enriched 
material has had only ~2 Myr to be expelled, mixed and 
diluted in the interstellar medium across these regions. 

In addition to the timing of these observations, we 
should also consider that our limited spatial resolution may 
be a factor. Even for nearby 30 Doradus that contains Of- 
and WR-type stars, nitrogen self-enrichment in that nebula 



has not been detected with any certainty (Rosa fe Mathis 



19871, and the H II region overall shows a normal N/O ra- 



tio for its 0.5 Zq metallicity ( Tsamis et al.|2003 
Pequignot|2005| . 



Tsamis & 



7 SUMMARY AND CONCLUSIONS 

We have analysed FLAMES-IFU integral field spectroscopy 
of Haro 11, a BCG well known for its high star-formation 
efficiencies and Lycv emission, making it analogous to the 
high-redshift LBGs. We have studied its morphology by cre- 
ating monochromatic emission line maps. The galaxy is seen 
to consist of three star-forming knots (named A, B and C), 
arranged in the shape of a heart with Knot B, the region 
of strongest Ha surface brightness in the centre and Knots 
A and C to the east and south-east, respectively. Haro 11 
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Knot A Knot B 

F{\) /(A) F(\) /(A) 



He I A4471 3.22±0.12 3.23±0.22 3.93±0.21 4.16±0.29 

He I A6678 3.10±0.21 3.09±0.23 3.55±0.10 2.89±0.11 

c(H/8) 0.01 ± 0.01 0.27 ± 0.02 

He+/H+ x 10 2 (A4471) 6.53 ± 0.41 8.31 ± 0.58 

He+/H+ x 10 2 (A6678) 8.34 ± 0.62 7.82 ± 0.28 

He+/H+ x 10 2 7.43 ± 0.54 7.74 ± 0.30 

ICF(He) 1.18 1.25 

He/HxlO 2 8.80 ± 0.64 9.68 ± 0.38 



Table 6. Regional He/H abundances for Knots A and B from integrated line fluxes (here in units relative to H/3 = 100) summed over 
the regions defined in Fig. [2] 






orcsecs 



0.00 0.30 0.59 0.89 1.18 1.48 
log(Flux x10" 



" i6 ergs s 1 cm" 2 ) 



Figure 11. Emission map of the blue WR feature, showing the 
distribution of WR stars throughout Haro 11. Overlaid dashed 
contours represent the Ho integrated flux. North is up and east 
is to the left. 




EW (H/3) / A 

Figure 12. Map of the equivalent width of H/3 across Haro 11. 
Overlaid red lines are contours of H/3-region continuum. North is 
up and east is to the left. 



exhibits complex emission line profiles, with most lines con- 
sisting of a narrow, central component, an underlying broad 
component, and a third narrow component. In agreement 
with previous studies of Ha kinematics, this system is not 
dynamically relaxed and is possibly an ongoing merger. In 
particular, the contrasting kinematics of the broad and nar- 
row emission line components suggest the presence of a 
counter-rotating disk within the central region of the galaxy. 

For the first time we present 200 pc resolution maps 
of T c , iV c , and the abundance ratios of O/H, N/H, Ne/H 
and S/H across Haro 11. The abundance maps were derived 
using the direct method of estimating electron temperature 
from the [O III] A4363/A4959 line ratio. 

Variations in oxygen abundance are seen across Haro 11 
in the narrow emission component, with a peak in abun- 
dance of 12 + log(0/H) = 8.25±0.15 in Knot B, decreasing 
to 8.09±0.23 and 7.80±0.13 in Knots A and C, respectively. 
We measured a difference of ~2000 K between the temper- 
ature of the broad and narrow line emission regions in Knot 
B, but the spectral resolution of our data is not sufficiently 
high to allow the decomposition of [O n] in the broad re- 
gions. Thus the total oxygen abundance differential between 
the two gas phases remains unknown. 

The low O/H abundance ratio in Knot C tallies with 
the increased T c in that region. We detect a high N/O ratio 



(previously reported for the whole galaxy Izotov & Thuan 
|1999| ) only in Knot C, with Alog(N/0)~ +0.5 dex, spatially 
coincident with the decrease in oxygen abundance. This knot 
also shows a distinctly lower N c in comparison to the rest of 
the galaxy. This may relate to the previously reported Lya 
emission observed solely from this region. 

We assess the reliability of abundances derived from 
galaxy-integrated spectra by comparing those derived from 
integrated spectra with spatial averages. Whilst abundances 
relative to oxygen are consistent within the uncertainties, 
abundances relative to hydrogen are not. These results illus- 
trate the superior information content of abundances from 
spatially resolved spectra. 

Using the 4686 A WR signature, we detect evidence of 
a large WR population in Knots A and B. No evidence for 
WR stars was found in Knot C, although we do measure 
an enhanced N/O and a decreased O/H ratio. This leads us 
to believe that the abundance anomaly in this region may 
be due to a recently completed WR phase where N remains 
undiluted by subsequent mass loss of lower mass stars. The 
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strong WR emission in Knots A+B combined with 'nor- 
mal' N/O levels and a young stellar population (<5 Myr) 
suggest we are observing these knots at an epoch where 
ejected N-rich material is yet to mix with the warm ionized 
gas. Alternatively, even the spatial resolution of these data 
(0.52" spaxeP 1 ), corresponding to ~210 pc at the distance 
of Haro 11), may be preventing us from detecting small-scale 
N-enhancements that have already formed. 

To conclude, the IFS data presented in this study show 
the complexity of Haro 11 as a system; not only is it com- 
plicated kinematically but also chemically. In order to bet- 
ter understand the evolution of this system, the inhomoge- 
neous conditions and metallicity of its ISM exposed by our 
chemodynamical study need to be combined with its star- 
formation history in models. Observations at even higher 
resolution are required to constrain the mixing timescales of 
the ejecta from its large WR population. Successful model- 
ing of an LBG analogue such as this could have a significant 
impact on our understanding of the chemical homogeneity, 
ISM mixing, and effect of feedback on galaxy evolution in 
both the local and high- z universe. 
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APPENDIX A: REGIONAL OBSERVED 
FLUXES AND LINE INTENSITIES 

Tables DTT1 1 A3I list the observed fluxes and de-reddened line 
intensities derived from summed spectra over the regions de- 
fined in Fig. [2] Fluxes are given for each separate velocity 
component (C1-C3), along with its de-reddened line inten- 
sity and FWHM. This paper has been typeset from a Tp^X/ 

DTpjX file prepared by the author. 
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23.41±1.01 


74.61±1.09 


28.68±0.52 


28.146±0.60 


259 


.11±2.11 


He I A6678 


2.85±0.31 


2.66±0.29 


67.35±4.69 


3.00±0.27 


2.95±0.27 


259 


.27±19.96 


[S n] A6716 


10.84±0.61 


10.13±0.59 


72.78±2.19 


15.37±0.35 


15.08±0.38 


239 


,06±2.35 


c(H/3) 




0.09 ± 0.01 






0.02 ± 0.01 






F(H/3)xl0 14 




1.99 ± 0.08 






2.73 ± 0.04 







C3 





F(A) 




7(A) 


FWHM 


[O n] A3727 










[O II] A3729 










[Nc in] A3868 










H8+He 1 A3889 










He 1 A4026 










[S 11] A4068 










H5 










H7 


43.92±8.93 


44. 


61±9.97 


58.48±5.17 


[O in] A4363 










He 1 A4471 










[Fc in] A4657 










H/3 


100.00±12.14 


100. 


,00±14.54 


57.76±2.51 


[O in] A4959 


171.58±18.76 


165. 


52±20.61 


59.28±2.15 


[O 1] A6364 










[N n] A6548 


5.57±0.71 


5. 


.35±0.75 


49.33±3.28 


Ha 


288.47±27.36 


276. 


.75±30.48 


62.37±1.11 


[N n] A6584 


12.48±1.30 


11. 


96±1.41 


42.66±1.71 


He 1 A6678 


3.78±0.76 


3 


,62±0.75 


58.65±5.99 


[S 11] A6716 


13.14±1.49 


13. 


14±1.84 


52.72±2.07 


c(H/3) 




0.06 


± 0.04 




F(H/3)xl0 14 erg s" 1 cm- 2 




0.86 


± 0.07 





Table Al. Haro 11 regional fluxes and dc-reddenened line intensities (both relative to H/3 = 100) and FWHMs (in krns" 1 ) for summed 
spectra over Knot A (as defined in Fig. [2J . Line fluxes were extinction-corrected using the c(H/3) values shown at the bottom of the table, 
calculated from the relative Ha, H/3 and H7 fluxes. 
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F(A) 




7(A) 


FWHM 




F(A) 






/(A) 




FWHM 


[O n] A3727 


;->44 


56 + 8 Q9 


4(;<) 


13+44.25 


256 


33±17.33 
















[O n] A3729 


355 


21±9.19 


483 


63±45.61 


256 


11±17.32 
















[Ne in] A3868 


87 


79±2.96 


115 


98±11.03 


212 


82±4.82 
















H8+He I A3889 


59. 


32±2.76 


78 


09±7.82 


209 


83±8.22 
















He I A4026 


4 


01±0.20 


5. 


09±0.51 


172 


,42±7.98 
















[O 11J 


4 


58 + 94 


5 


754-0 58 


188 


22±8.86 
















H<5 


24 


43±1.31 


30 


42±3.07 


121 


96±2.80 


24 


,06±0 


.60 


24 


76±1.72 


296 


.04±4.70 


Fc II A4286 


1 


69±0.32 


2 


00±0.41 


213 


13±42.98 
















H 7 


39 


50±2.51 


46 


04±4.75 


121 


30±3.00 


46 


.56±1 


.23 


47 


50±3.19 


298 


,60±4.44 


[O in] A4363 


3. 


60±0.66 


4 


17±0.83 


132 


66±17.21 


2 


69±0 


.57 


2. 


68±0.57 


412 


.75±110.16 


Ho I A4471 


4. 


78±0.46 


5 


37±0.67 


122 


62±6.51 


3 


,51±0 


.22 


3. 


56±0.31 


384 


.48±28.20 


[Fc III] A4657 


4. 


47±0.28 


4 


75±0.46 


204 


16±11.96 
















H/3 


100. 


00±3.57 


100 


00±8.06 


118 


66±1.63 


100 


,00±1 


.63 


100 


00±5.71 


290 


.99±2.21 


[O in] A4959 


118 


11±3.62 


114 


74±8.82 


104 


61±1.16 


116 


,51±6 


.43 


111 


89±4.68 


276 


68±1.47 


[O I] A6364 


2. 


06±0.20 


1 


47±0.16 


119 


,18±5.77 


1 


99±0 


10 


1. 


91±0.12 


325 


.55±16.60 


[N II] A6548 


13 


61±0.61 


9. 


28±0.62 


89 


60±1.55 


14 


,41±0 


.28 


13. 


71±0.58 


200 


,66±2.06 


Ho 


417 


83±10.80 


286 


85±16.06 


125 


02±0.41 


306 


,78±3 


.67 


292 


02±11.54 


294 


,78±0.41 


[N II] A6584 


56. 


97±1.59 


38. 


59±2.17 


100 


,03±0.70 


55 


,54±0 


.69 


52. 


77±2.06 


264 


.94±1.21 


He I A6678 


2. 


80±0.22 


1 


87±0.17 


91 


73±3.41 


3 


,92±0 


10 


3. 


72±0.16 


242 


,46±6.23 


[S II] A6716 


15 


05±0.50 


9 


96±0.58 


97 


79±1.26 


25 


,50±0 


.33 


24 


16±0.93 


253 


,33±1.00 


c(H/3) 






0.52 ± 0.03 












0.07 ± 0.02 






F(H/3)xl0 14 erg s _1 cm~ 2 






5.83 ± 0.15 












11.88 ± 0.14 







Table A2. Same as for Table [Al] for summed spectra over Knot B (as defined in Fig.[2|. 



Knot C CI C2 

F(A) 1(A) FWHM F(A) /(A) FWHM 



[O II] A3727 


208. 


00±9.85 


245 


20±28.29 


125 


59±0.72 












[O II] A3729 


263. 


87±12.46 


311. 


07±35.88 


125 


47±0.72 












[No III] A3868 


53. 


77±4.66 


62 


38±8.40 


125 


33±10.34 












H8+He I A3889 


19 


62±4.75 


22 


71±5.98 


81 


12±21.74 












Ho I A4026 


1 


20±0.40 


1. 


36±0.47 


114 


56±40.39 












[S II] A4068 


6. 


99±0.63 


7. 


89±1.06 


135 


69±11.52 












US 


37. 


01±1.84 


41 


60±4.61 


99 


87±1.73 












Fo II A4244 


6. 


65±1.05 


7 


32±1.36 


174 


51±26.47 












Fc II A4286 


8. 


56±1.21 


9. 


36±1.59 


188 


33±26.61 












H 7 


43 


32±7.81 


47 


01±9.57 


80 


67±5.68 


43 


84±9.31 


47.54±11.10 


172 


,37±21.11 


[O in] A4363 


6. 


77±1.98 


7. 


32±2.24 


185 


05±59.63 












He I A4471 


3. 


86±0.39 


4 


10±0.56 


82 


76±8.15 












[Fe III] A4657 


5. 


37±0.75 


5 


55±0.92 


185 


77±26.72 












H/3 


100. 


00±6.58 


100 


00±10.65 


82 


08±2.13 


100. 


00±9.40 


100.00±12.74 


200 


,87±6.83 


[O in] A4959 


109. 


67±6.29 


107 


99±10.79 


81 


02±1.60 


101 


58±8.18 


100.04±11.63 


213 


,87±6.95 


[O I] A6364 


3. 


08±0.23 


2. 


58±0.25 


111 


20±6.54 












[N II] A6548 


11 


81±1.39 


9 


62±1.26 


67 


18±3.16 


17 


62±1.95 


14.39±1.80 


147 


,45±7.96 


Ha 


355. 


13±16.74 


290 


59±21.64 


97 


37±0.58 


357. 


83±23.94 


293.34±26.19 


286 


,36±1.29 


[N II] A6584 


50. 


24±2.59 


40 


82±3.13 


78 


31±1.09 


50. 


36±3.58 


41.00±3.77 


202 


,66±4.16 


He I A6678 


4. 


09±0.24 


3 


29±0.26 


88 


60±3.08 












[S II] A6716 


28. 


18±1.60 


22 


62±1.79 


79 


02±1.38 


23 


39±1.94 


18.81±1.89 


191 


,49±5.89 


c(H/3) 






0.28 ± 0.04 










0.28 ± 0.04 






F(H/3)xl0 14 erg cm -2 






1.07 ± 0.05 










0.81 ± 0.05 







Table A3. Same as for Table [Al] for summed spectra over Knot C (as defined in Fig. [2|. 



